Efficient solvability of Hamiltonians and limits on the power of some quantum computational 

models 



Rolando Somma, Howard Barnum, and Gerardo Ortiz 
Los Alamos National Laboratory, MS D454, Los Alamos, NM 8754^ 

Emanuel Knill 

Mathematical and Computational Sciences Division, National Institute of Standards and Technology, Boulder CO 8030^ 

(Dated: February 1, 2008) 

We consider quantum computational models defined via a Lie-algebraic theory. In these models, specified 
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acts. Similar results hold for the computation of the expectation value of operators implemented by a gate- 
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Quantum models of computation are widely believed to be 
more powerful than classical ones. Although this has been 
shown to be true in a few cases, it is still important to deter- 
mine when a quantum algorithm for a given problem is more 
resource efficient than any classical one, or, conversely, when 
a classical algorithm is just as efficient as any quantum coun- 
terpart. In general, one needs to know whether it is worth 
investing in building a quantum computer (QC) and what is 
required for success. In this paper, we show close connec- 
tions between these issues and the efficient (or exact) solv- 
ability of Hamiltonians. In particular, we show that a class of 
quantum models we call generalized mean-field Hamiltoni- 
ans (GMFHs) iQl is efficiently solvable and furthermore does 
not provide a stronger-than-classical model of computation: 
A quantum device engineered to have dynamical gates gener- 
ated by Hamiltonians from such a set cannot directly simulate 
universal efficient quantum computation and can be efficiently 
simulated by a classical computer (CC). 

An algorithm is a sequence of elementary instructions that 
solves instances of a problem. It is said to be efficient if the 
resources required to solve problem instances of size N are 
polynomial in N (poly(iV)) resources. Typically, the size of 
a problem instance is the number of bits required to represent 
it, and the relevant resources are time and space. In the last 
few years it has been shown that many pure-state quantum al- 
gorithms can be efficiently simulated on a CC when the extent 
of entanglement is limited (e.g., j^lH) or when the quantum 
gates available are far from allowing us to build a set of uni- 
versal gates |i4||5||6|i|. Here, we focus on a Lie algebraic anal- 
ysis to obtain other situations where quantum algorithms can 
be efficiently simulated by CCs. The so-called generalized 
coherent states (GCSs) 01 play a decisive role in our analysis. 

The algorithms considered here make use of the Lie- 
algebraic model of quantum computing (LQC). An LQC algo- 
rithm begins with the specification of a semisimple, compact 



7\ /-dimensional real Lie algebra f) of skew-Hermitian opera- 
tors acting on a finite-dimensional Hilbert space H, with Lie 
bracket [X, Y] := XY - YX. Without loss of generality, the 
action is irreducible. The algorithm begins with a maximum- 
weight state |hw) in H and applies gates expressed as expo- 
nentials for certain X int). The output of the algorithm is 
a noisy expectation of an operator in f) or in e''. LQC algo- 
rithms cannot trivially be classically simulated because of the 
possibility that the dimension of Ti is exponential in the speci- 
fication complexity of f) and | hw) . In order to precisely define 
the model of LQC we require some results from the theory of 
Lie algebras. See Q for a textbook covering the basic theory 
of Lie algebras. 

Our intention is to restrict observables and Hamiltonians to 
operators in -y/— If) . The dimension of H may be exponential 
in AI. Since we wish to implement computations with re- 
sources that are polynomial in M, our knowledge of t) cannot 
involve explicit matrix representations of its operators. We 
therefore assume that f) is specified as an abstract Lie alge- 
bra [) together with a "maximum weight" w characterizing 
its action on H. For computational purposes, we also use a 
small-dimensional faithful representation of i). To be specific, 
we use the adjoint representation, but for efficiency, one can 
choose the first fundamental representation instead. We use 
the following notational conventions; Objects with a "hat" 
C) belong to the representation of t) on H. Objects with an 
"overline" C) belong to the chosen faithful representation. Lie 
algebraic objects with neither a hat nor an overline are associ- 
ated with the abstract Lie algebra (representation unspecified). 
Implicit in these conventions are the representational isomor- 
phisms i) ^ i) and f) f). 

For the purpose of efficient representation, it is convenient 
to work with the complexification Cf) of f) and use a Cartan- 
Weyl (CW) basis (see, for example, 0]) for Cf). Thus, we 
assume a decomposition Cf) = i) d ® (B i)^ , where f|_D is 
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a Cartan subalgebra (CSA), and [)+ and t)~ are algebras of 
generalized raising and lowering operators, respectively, f) d 
is linearly spanned by named elements hi, . . .h^, and f}^ by 
e^_^ , . . . , e . The aj are linear functionals on ijo called the 
positive roots of The abstract Lie algebra is specified by 
the identities [hk,hi] = 0, [/ife,e± ] = ioj-feC^ , [e+.e"] = 
J2k hjhk, and for j ^ k, [e+. , e±J = Cjfce+.±„^ . The bases 
of and ()^ may be chosen so that the "structure constants" 
ttkj, bkj and Cjk are ratios of integers with poly(i\/) digits. 
The structure constants do not uniquely specify the action of 
t) on H. According to the representation theory of semisim- 
ple complex Lie algebras, this action is uniquely specified by 
its "maximum weight", which is a Unear functional w on [)£) 
given by its values w{hk) on the distinguished basis of \)d- 
The w{hk) are integral and are the eigenvalues of hk on the 
unique state |hw) annihilated by [)+: /ife|hw) = w{hk)\hvj)- 
The state |hw) is called the "maximum weight state" of the 
representation and its orbit under e'^ generates the family of 
GCS of Ti, with respect to (). The Hermitian inner product of 
Ti. and the Hermitian transpose of operators on Ti. induce a cor- 
responding Hermitian transpose operation on Cf]. We assume 
that the CW basis is chosen so that the Hermitian transpose is 
given by h\ = hk and {e^. )t = e^^. . We also assume that the 
linear space on which f) acts is endowed with a Hermitian in- 
ner product for which the representation [) is skew-Hermitian 
and the Hermitian transpose matches the one defined for f). 

The formal specification of an LQC algorithm requires the 
structure constants of an abstract CW decomposition of [) and 
the weight coefficients w{hk) determining |hw). The spec- 
ification complexities of f) and | hw) are the number of bits 
required to represent the numerators and denominators of the 
structure constants and the w{hk)- Thus they are polynomial 
in the dimension of f) and logmax(w(/ife)). The gates of the 
algorithm may be unitary exponentials e*'^ , with X a mem- 
ber of the CW basis. The gate's resource requirement is the 
number of bits required to represent t plus \t\. More gener- 
ally, we can allow as gates any with H E ij, where the 
resource requirement is given by the specification complexity 
of (defined below). There are several alternatives for how 
the algorithm's output is obtained. We consider two. In the 
first, the output is obtained by measuring the expectation of 
an operator A G Cf). In the second, it is obtained from the 
absolute value of the expectation of an operator U G e''''. The 
resource cost of making the measurement is proportional to 
the sum of number of bits of precision and the specification 
complexity of A or U. The specification complexity of A is 
that of A (the corresponding operator in the abstract Lie al- 
gebra f)) and is given by the number of bits used to represent 
the coefficients of A when expressed in the CW basis. If U is 
of the form with H E Cf), its specification complexity is 
that of H plus maxdi/al) where Ha ranges over the coeffi- 
cients of H expressed in the CW basis. Our assumption about 
the resource cost of measurement (i.e., the number of bits b 
of precision) makes the LQC model just defined very power- 



ful but physically unreasonable. In particular, an LQC algo- 
rithm gives exponentially better precision than an algorithm 
of similar resource cost for the standard quantum computa- 
tional model. The standard quantum algorithm would need to 
be repeated exponentially many times in b to return an expec- 
tation value with b bits of precision. In the standard quantum 
computational model, the hypothetical ability to determine ex- 
pectation values with b bits of precision using resources poly- 
nomial in b implies the ability to efficiently solve problems 
in i^P, the class of problems associated with the ability to 
count the number of solutions to NP-complete problems such 
as satisfiability. This is a consequence of more general results 
in 

A natural question is when and how LQC can efficiently 
simulate, or be simulated by, standard quantum or classical 
computation. The measurement models we introduced for 
LQC have the same form as many typical problems in physics, 
which involve the evaluation of correlation functions 

{W)^Tr[pW], (1) 

where p ~ J2s=i Pali's) {4's\ is the density operator of the 
system (ps > 0',J2sPs = 1)' \4's) are pure states, and W is 
a Hermitian or unitary operator acting onH. In general, the 
dimension dofH increases exponentially in the problem size 
N , where the problem size is determined by quantities such as 
the volume or number of particles of the system. An algorithm 
to evaluate (W) with accuracy e is efficient if the amount of 
resources required is bounded by polylog(d) + poly (1/e). 

An efficient quantum algorithm to evaluate Eq. Q exists 
if the state p (or a good approximation to it) can be effi- 
ciently prepared on a QC and if W can be efficiently mea- 
sured by usi ng, f or example, the indirect techniques described 
in Refs. HUfiai . Unfortunately, known classical algorithms 
for this purpose typically require resources polynomial in the 
dimension d, which can be exponential in the problem size N . 
However, if the problem can be specified Lie algebraically, 
this classical complexity can be greatly reduced and exponen- 
tial rather than polynomial accuracy is efficiently achievable. 

Theorem 1 With p as defined following Eq. 0, if \4)s) = 
e^=|hw) are GCS's oft) (A, € I)) and W € Ct), then (W) 
can be classically computed to accuracy e in time polynomial 
in log{\/e) and the sum of the specification complexities ofi), 
|hw), W, As and ps. 

Proof: We have (W) = ^^^-^Ps(hw|Ws|hw), where 
Ws = e-^^We^^. In the CW basis, 

r I 

ws^Y. + E + ^r««. ' (2) 

k=l 3 = 1 

where u^, f G C. To obtain these coefficients, we can com- 
pute Ws in the adjoint representation: Ws = e'^'We^" = 

r _ I 

J2 u^hk + X) + v^~e~.. To compute the ul and 

k=l 3 = 1 
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Vj to accuracy 6 requires computing the matrix exponentials 

e^^' , and matrix multiplication followed by an expansion of 
the resulting matrix in terms of the CW basis. The matrix ex- 
ponentials can be obtained to accuracy S' (in the 2-norm) in 
time polynomial in log(l/5') and the maximum of the entries 
of the Ashy direct series expansion or other, more efficient 
methods 11311 . Matrix multiplication and basis expansion in- 
crease the 2-norm error by at most a constant factor, so that the 
ul and Vj can be efficiently obtained to the desired accuracy. 

Using the property that the e^^ either map |hw) to an or- 
thogonal state or annihilate it, we rewrite Eq. Q as 

L r 

(W^) (3) 

s=l k=l 

and this sum can be evaluated efficiently with respect to the 
given specification complexities. ■ 

The following variant of Thm.^holds for W — with 
H e C[). 

Theorem 2 = e^=|hw) (A^ e i)) are GCS's oft) and 

W = with H £ Cf), then \ {W)\'^ can be classically com- 
puted to accuracy e in time polynomial in log{l/e) and the 
sum of the specification complexities of l), |hw), W, As and 

Ps- 

Proof: We can expand | (W) p as 

^ y^^PsPs'tr ds,s', with 6s, s' := 

s.s' 

|hw)(hw|e-^=e-^e^=|hw)(hw|e-'^^'e^^^='|hw)(hw| . 

Os.s' is proportional to |hw) (hw| and its trace is the constant 
of proportionality. We can express |hw)(hw| as a limit of 
operators in e^''. Let L = X]I-=i w{hk)hk and define oj by 
i|hw) = tj|hw). Then {il}\L\il}) < uj for ^ |hw), from 
which it follows that |hw)(hw| lim^^oo e"*"e*^ Q- Be- 
cause the eigenvalues of L are integral, convergence is expo- 
nentially fast in t. Let 

E(t) - ^p,p,,e-3-*e*^e-^=e^e^=e*^e-^='e^'e^^'e*^. 

E(t) is positive definite Hermitian and converges to Os.s' as 
t ^ oo. For a given t, we can compute E{t) by computing ex- 
ponentials and multiplying matrices in the adjoint representa- 
tion. Observe that the maximum eigenvalue nit) of E(t) con- 
verges exponentially fast to | (yV) p. To compute nit) we first 
determine Q(t) such that Eii) = e'^^*). With the assumed 
Hermitian inner product on the adjoint representation, E(t) is 
positive definite. Thus, there is a unique Hermitian Q(t) sat- 
isfying E{i) = e'^^*), and Q{i) is necessarily in \J —Vc). The 
operator Q{t) can be obtained via any conventional efficient 



diagonalization procedure for non-negative definite matrices. 
We can then use an efficient Jacobi-like diagonalization proce- 
dure fl^ to obtain unitary operators Uit) € ^ and q(t) Ghjj 
such that C7(t)g(<)C7(<)t = Q{t). The maximum eigenvalue 
of E{t) is given by the exponential of the maximum eigen- 
value of q{t). At this point we require a number of results 
from the representation theory of Lie algebras. For example, 
see II2J. The element q{t) induces an alternative order on the 
roots, according to which a root aj is positive if aj{q{t)) is 
positive. (To remove degeneracies, it may be necessary to 
slightly perturb q{t).) For this ordering, we determine sim- 
ple roots Pk and corresponding members h'j^ G t)D such that 
/I'j, is isomorphic to hk via a member of the Weyl group. We 
can expand q{t) = qjh'j. Uniqueness of maximal weights 
in representations of Lie algebras implies that the maximum 
eigenvalue of q{t) is given by w'{q{t)) = J2j w{hk)qj- 

We claim that the necessary steps can be implemented with 
polynomial resources in the dimension of the Lie algebra and 
the number of digits of precision of K{t). The matrix and 
root manipulations can be implemented efficiently, but with 
respect to the precision of entries of the matrix. It is neces- 
sary to realize that unless the weight w is sufficiently small, 
E{t) converges to exponentially fast in t. However, be- 
cause the are polynomial in M, the rate of conver- 
gence to is bounded by e^P°'^(*'^^ To compute n{t) to a 
desired number P of digits of precision, it suffices to compute 
in the low dimensional matrix representation with a precision 
of poly(A/) + poly(P) digits, which can still be done with 
polynomial resources. The relevant Weyl group transforma- 
tions can be done efficiently by use of one of the constructive 
proofs of the transitivity of the Weyl group. See, for example, 
the proof of Thm. 2.63 in j^l- ■ 

Important special cases motivating these results are 
fermionic linear optics quantum computation (and equivalent 
matchgate models introduced by Valiant), which is efficiently 
classically simulatable 11811 . and models that also in- 

clude linear fermionic operators (so(2A^ + 1)) for which an 
extension of the canonical Bogoliubov mapping exists ITgll . 
Natural bosonic analogues of the fermionic results also exist. 
As a result it is possible to efficiently simulate quantum com- 
putational models in which coherent states are acted on by 
linear optical circuits, and measured via homodyne detection 
0], and of models with initial multimode squeezed states and 
squeezing gates as well as hnear ones I20I1 . Like LQC with the 
second measurement strategy, these involve the efficient sim- 
ulation, in the dimension of a Lie algebra, of a computational 
model in which coherent states of a Lie group with gates gen- 
erated by the algebra constitute the initial states and computa- 
tion. However, in the bosonic case the relevant algebra is not 
semisimple, and the relevant irreps are infinite-dimensional. 

We can now address the important question of the classical 
simulatabihty of LQCs. 

Theorem 3 For both LQC measurement schemes, the result 
of an LQC algorithm A can be obtained by use of classical 
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computation in time polynomial in the specification complex- 
ity of A 

Proof: The action of the gates of the algorithm result in the 
state 10) = nm=i ^'^'"Ihw) where e ^. Let 

t t 
(Vi^) = (hw| e-^"VF e^-lhw). (5) 

m— 1 m—1 

The result of the algorithm is {W) if e t). or | {W) \ if 
W G e"''' . The result can be computed by generalizing the 
algorithms given in the proofs of Thms. ^ and |2 All that 
is required is to compute the full product nm=i ^^"^ instead 
of the single exponential required for Thms. ^ and |2] The 
complexity of the method is polynomial with respect to the 
specification complexity of A. H 

The meaning of Thrnjsjcan be expressed in terms of gen- 
eralized entanglement Illll4i EUl . A pure state is generalized 
unentangled (GU) with respect to a preferred set O of observ- 
ables if it is extremal among states considered as linear func- 
tionals on O, otherwise it is generalized entangled. In a Lie 
algebraic framework, a GU state is a GCS of a semisimple 
compact Lie algebra. Thus, Thm. |3] states that if a quantum 
computation does not create generalized entanglement with 
respect to a polynomial-dimensional semisimple compact Lie 
algebra, such a computation can be efficiently simulated to 
exponential precision on a CC. We also remark that, because 
it cannot access all pure states during the computation, but 
only the submanifold of generalized unentangled ones, such a 
computation cannot directly simulate standard quantum com- 
putation. 

For appUcations to physics simulation, the following corol- 
lary shows that higher-order correlation functions can also be 
computed efficiently, provided the order is not too large. 

Corollary 1 Let W^, . . . , W'^ be operators in <C\). For fixed 
q, the expectation value of correlation functions of the form 
(W^-'-VV) = J2s=iPs{^Me~^'W^ ■ ■ - Wie^^lhw), can 
be computed on a CC in time polynomial in log(l/e), and the 
sum of the specification complexities oft), |hw), W\ As and 

Ps- 

The complexity of our algorithm for computing the corre- 
lation function in the corollary is exponential in q. 

Proof: We outline an efficient algorithm for computing the 
desired correlation function. First we expand each VKf = 
g-Asyyj j-jjg basis as in the proof of Thm. ^ The 
desired correlation is given by X]s(^*l Yij W^\hw). We for- 
mally multiply the CW basis expressions for the to obtain 
sums Ps of formal products of members of CW basis rep- 
resenting the Ylj W^. Each product of CW basis members 
is standardized by using the commutation rules so that each 
term is a product where all lowering (raising) operators follow 
(precede) members of This is similar to the procedure of 
Wick's theorem. After this transformation, terms that retain 



some lowering or raising operators contribute nothing to the 
correlation functions. The remaining terms' contribution is 
easily computed from (hw|/ife|hw) = w{hk). The contribu- 
tion to the complexity of the procedure of the formal multipli- 
cation and standardization procedure grows exponentially in 
q. The number of terms that arise is bounded by poly(Af)'', 
so that for fixed q, the complexity remains polynomial in the 
given specification complexities. Further details are available 
in m. ■ 

The algorithms given above can also be used to analyze 
certain interacting physical models. We use the term GMFH 
fll for Hamiltonians belonging to V~lf) for f) in a sequence 
of semisimple compact operator Lie algebras of dimension 
M < polylog((i) acting on d-dimensional Hilbert spaces. A 
GMFH is necessarily specified in terms of a basis of f) that 
can be efficiently transformed to a CW basis. An example 
of a GMFH is given by the N spin- 1/2 Ising model in a 

transverse magnetic field Hj = ^ {galai^^ + ""I), where 

Hi is an element of the Lie algebra 5o(2A^), with dimension 
M = 2N^ - N = polylog(rf), where d = 2^. Interestingly, 
this model can be exactly solved and, as we will show, this 
result can be extended to any GMFH. We say that a Hamilto- 
nian acting on a d dimensional Hilbert space can be efficiently 
(exactly) solved when any one of its eigenvalues and a de- 
scription of the corresponding eigenstate can be obtained and 
represented to precision e in polylog((i) + poly(l/e) compu- 
tational operations on a CC. In general, this definition makes 
sense when we focus on Hamiltonians describing the inter- 
actions of iV-body systems, where d increases exponentially 
with N. 

Theorem 4 GMFHs can be efficiently solved. 

Proof: Let Hi\if be a GMFH in \/— Tf) given in terms of 
a CW basis of (j as in Eq.|2l We show that to solve Hmf it 
suffices to diagonalize it according to 

r 

Hd = UHmfU^ = Y.^khk, (6) 
fe=i 

with ffe G R and [/ G e'' unitary. The eigenvalues of H]\if are 
shared with those of H]j. A description of the corresponding 
eigenspaces consists of an eigenspace of transformed by 
U^, where U may be described by a sequence of LQC gates. 
According to the representation theory of Lie algebras, the 
eigenspaces of Hd consist of weight states of t), which can 
be obtained from the highest weight state by applying lower- 
ing operators. They are characterized by linear functionals A 
on Hu of the form X{hk) = w{hk) — X); '^i(^i{hk), where 
the 71; are non-negative integers. Which choices of ?i; cor- 
respond to weight states is readily determined from the rep- 
resentation theory of Lie algebras. Once we have expanded 
H]j ~ J2k ^khk, the eigenvalue corresponding to A is readily 
computed as X{Hd) = J2k SkHhk). 
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To efficiently diagonalize Hmf and obtain a specification 
of U, we compute in the adjoint representation and apply a 
generalization of the Jacobi method 115112311 to Hmf- It yields 
an exponentially converging diagonalization and an expres- 
sion for U in terms of a sequence of exponentials of members 
of the su(2) subalgebras of Ct) generated by the pairs e^.. 
This suffices for our purposes. ■ 



Example. 

Si, 7=1 *y ("^iC 



The fermionic 



Hamiltonians Hm f = 
h.c, where the oper- 



ator c] (Cj) creates (annihilates) a spinless fermion at the ith 
site, belong to a representation of the Lie algebra so{2N) of 
dimension M = 2N'^ — N < polylog((i). A faithful represen- 



, and c. c 



Tn+i 



-N+j,N+i, 

where 



have H-1 in the kth row and fc'th 



tation of so(2iV) is given by c; c 

the 27V X 2N matrices Tkk 
column, and zeros otherwise. Therefore, we write the matrix 
of Hm f in this representation and apply the Jacobi algorithm 
to diagonalize it. The result is equivalent to the one given by 
the Bogoliubov transformation 12411 . where the Hamiltonian 
maps as Hmf ^ Hd = J^l^i ^k{lllk ~ 1/2), where the 
operator (7^) creates (annihilates) a fermionic quasiparticle 
in the fcth mode. 

Although LQC algorithms and GMFHs can be efficiently 
simulated or solved on a CC, it may still be useful to imple- 
ment the algorithms or simulate GMFHs with QCs. In par- 
ticular, there may be problems where a key component is ex- 
pressed in terms of LQC or GMFHs but a more complex quan- 
tum computation is required to determine the information of 
interest. One case of interest is where the LQC or GMFH 
component requires preparing a GCS. One way for such a 
GCS to arise is as the ground state of a GMFH. According 
to the next theorem, such GCSs are efficiently preparable on a 
QC that has efficient access to the LQC initial state and gates. 

Theorems Let a GCS 10) 0/ ft be specified as the ground 
state of a Hamiltonian H G Then \4>) can be prepared 

by use of resources polynomial in the specification complexity 
of H on a QC with the ability to initialize |hw) and efficiently 
apply LQC gates. 

Proof: It suffices to determine a. U £ expressed as a 
polynomial product of LQC gates such that H = UHuU^ 
with E l)D such that D induces the root order associ- 
ated with the CW basis. (See the proof of Thm.|2]for how an 
element of 1)d induces a root order.) The state |0) is then ob- 
tained as J7| hw) and hence is efficiently preparable using LQC 
operations. To determine U we can first use the generalization 
of the Jacobi method as discussed previously. This yields an 
element of f)^) that does not necessarily induce the desired 
root order To complete the determination of U requires using 
a sequence of Weyl reflections to obtain the desired root order 
The sequence may be obtained using the method mentioned at 
the end of the proof of Thm|2] ■ 

Our results provide an alog ues of the Gottesman-Knill the- 
orem 101 (cf. also I25I 13) concerning the efficient sim- 



ulatability of Clifford-group computational models, and of 
results on the simulatability of certain multimode coherent- 
state and squeezed-state computational models 15 Eoll . One 
might hope for a treatment, perhaps based on Lie groups and 
groups of Lie type, that will unify these results, specifically 
those based on (1) finite dimensional semisimple Lie alge- 
bras, (2) Bosonic linear optics with homodyne detection (tied 
to an infinite-dimensional irreducible representation of a solv- 
able Lie algebra) and possibly squeezing (involving a nilpo- 
tent Lie algebra), and (3) Clifford groups and semigroups. 
Our results cast additional light on why quantum computers 
may be more powerful than classical computers. It is a cru- 
cial fact that the generators of its gate-set, though their num- 
ber can be chosen to grow polynomially, generate an expo- 
nentially large Lie algebra acting on an exponentially large 
Hilbert space. If the growth of the dimension of the gener- 
ated Lie algebra is polynomial, a computation with this gate 
set using compatible state preparations and measurements can 
be simulated with polynomial efficiency on a classical com- 
puter by working in a low-dimensional faithful representation 
of the Lie algebra. What other algebraically constrained mod- 
els of quantum computation are efficiently classically simu- 
latable? Such structures may underlie the efficient solvability 
of further classes of Hamiltonians of condensed matter mod- 
els, which go beyond the GMFHs, such as those solvable via 
a Bethe-type Ansatz. 
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Ref. i;15] and the US DOE and NSA for support. Contribu- 
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